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DERIVATION OF THE NLS BREATHER 
SOLUTIONS USING DISPLACED 
PHASE- AMPLITUDE VARIABLES 

Natanael Karjanto and E. van Groesen 

Abstract. Breather solutions of the nonlinear Schrodinger equation are derived in this 
paper: the Soliton on Finite Background, the Ma breather and the rational breather. A 
special Ansatz of a displaced phase-amplitude equation with respect to a background is 
used as has been proposed by van Groesen et. al. (2006). Requiring the displaced phase 
to be temporally independent, has as consequence that the dynamics at each position 
is described by the motion of a nonlinear autonomous oscillator in a potential energy 
that depends on the phase and on the spatial phase change. The relation among the 
breather solutions is confirmed by explicit expressions, and illustrated with the amplitude 
amplification factor. Additionally, the corresponding physical wave field is also studied 
and wavefront dislocation together with phase singularity at vanishing amplitude are 
observed in all three cases. 

Key words and Phrases: displaced phase-amplitude variables, nonlinear Schrodinger equation, 
breather solutions, amplitude amplification factor, wavefront dislocation. 



1. INTRODUCTION 

The nonlinear Schrodinger (NLS) equation is a nonlinear dispersive partial 
differential equation that has been used as a mathematical model in many areas of 
Mathematical Physics, such as nonlinear water waves, nonlinear optics and plasma 
physics (Ablowitz and Segur, 1981; Ablowitz and Clarkson, 1991). In this letter, 
we will consider the spatial NLS equation, given as follows: 

+ ifldfy + z 7 |^| V = 0. (1) 

We will derive and study three exact solutions of this equation. These exact solu- 
tions are known in the literature as 'breather' solutions of the NLS equation (Dysthe 
and Trulsen, 1999; Dysthe, 2001; Grimshaw et. al., 2001). The name 'breather' 
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reflects the behavior of the profile which is periodic in time or space and localized 
in space or time. The concept was introduced by Ablowitz, Kaup, Newell and Se- 
gur (Ablowitz, et. al., 1974) in the context of the sine-Gordon partial differential 
equation. The breather solutions are also found in other equations, for instance in 
the Davey- Stewart son equation (Tajiri and Arai, 2000) and modified Korteweg-de 
Vries (mKdV) equation (Drazin and Johnson, 1989). 

To find the breather solutions of the NLS equation, we use a special Ansatz in 
the description with displaced phase-amplitude variables introduced by van Groesen 
et. al. (2006), with a displacement that depends on the background. Requiring 
the displaced phase to be temporally independent, has as consequence that the 
dynamics at each position is described by the motion of a nonlinear autonomous 
oscillator in a potential energy that depends on the phase and on the spatial phase 
change. It is remarkable that this assumption leads to the breather solutions of the 
NLS equation that are known in the literatures: the Soliton on Finite Background 
(SFB) (Akhmediev et. al., 1987), the Ma breather (Ma, 1979) and the rational 
breather (Peregrine, 1983). These solutions have been suggested as models for a 
class of freak wave events seen in 2 + 1 dimensional simulations of surface gravity 
waves (Henderson et. al., 1999). However, the SFB shows many properties of 
extreme, or freak, rogue wave events, as observed by Osborne et. al. (2000), Calini 
and Schober (2002). The SFB is also a good candidate in order to generate extreme 
wave events in a hydrodynamic laboratory (Andonowati et. al., 2006; Huijsmans 
et. al., 2005). 

This letter is organized as follows. Section 2 presents some exact solutions 
of the NLS equation. We demonstrate that three different solutions of the NLS 
equation are obtained within a similar formulation. Section 3 explains briefly the 
relationship among the solutions of the NLS equation. Section 4 will discuss the 
evolution of breather solutions in the Argand diagram. In Section 5 we will deal 
with surface water waves and observe the occurrence of wavefront dislocations in 
the corresponding density plots of the breather solutions. The final section draws 
some conclusions and remarks on the subject in this letter. 



The simplest nontrivial solution of the NLS equation is the plane-wave or 
the 'continuous wave' (cw) solution. It does not depend on the temporal variable 
r and is given by Ao(£) = roe~* 7r °^. Another simple solution is the 'one-soliton' 
or 'single-soliton' solution. It can be found by the inverse scattering technique 
(Zakharov and Shabat, 1972), but also much simpler by seeking a travelling- wave 
solution. An explicit expression is given by 
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It will become clear when we consider the corresponding physical wave field that 
both the cw and the one-soliton solutions have coherent structures, while other 
solutions that we will deal with below do not. 

In van Groesen et. al. (2006), a displaced phase-amplitude description has 
been introduced. As a special case, we will look for a solution of the NLS equation 
in the form 

A(t,T)=A (0[G(l;,T)e i ^ -l], (3) 

for which the phase </>, which is 'displaced' since the background contribution is 
taken out, is taken to depend only on the spatial variable £. Substituting ([3| into 
the NLS equation 0, we obtain two equations that correspond to the real and the 
imaginary parts, respectively. 

Multiplying the real part by cos(/> and the imaginary part by shi(/>, adding 
both equations will give a Riccati-like equation for G: 

d^G + 7^ sin 2(j)G - 7^ sin (j)G 2 = 0. (4) 

Solving this equation leads us to a conclusion that G can be written in a special 
form. By letting G = then equation Q becomes a first order linear differential 
equation in H: 

d^H - 77*0 sin 2<j) H + 7^ sin = 0. (5) 

Let us take P(£) = exp (—77*0 f sin20(£) d£) as the integrating factor, multiply it 
to ([5| and integrate the result with respect to £, then we obtain the solution for H: 

H(£ = -irl J P(Q sin m - Or) = Q(Q - C(r) 

where — ((r) is a constant of integration that depends on r and Q(£) denotes 
the term with the integral sign. Assuming that the displaced phase is an 
invertible function, we can write £ = £(0) and dropping the tilde signs to indicate 
that P = P(<j>) and Q = Q ((/)). Consequently, G is now written as follows: 

We will observe in the following subsections that by choosing three different func- 
tions of C(t) will lead to three different breather solutions of the NLS equation. 

On the other hand, multiplying the real part by sin (j) and the imaginary part 
by cos 0, subtracting one from the other we obtain a nonlinear oscillator equation 
for G: 

pd 2 T G + + 2 7 r^ cos 2 <j>)G - 3 7 r 2 cos 0G 2 + 7 r 2 1 G 3 = 0. (8) 

In deriving the breather solutions, we will compare this equation with the second 
order differential equation for G after choosing a particular function of C( r )- 
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Firstly, we take £(r) = cos(^r), where v is a modulation frequency. For the 
normalized quantity v = vj {^Q^J~^j we consider < v < y/2. The differential 



equation for G becomes: 

d 2 G = -v 2 G + iv 2 %G 2 + 2v 2l -^$^G\ 

P P Z 

By comparing with (J8|, the solution (JT]) is obtained with P(4>) = v 2 Q((j))/ cos0, 
Q 2 ((j)) = 2 cos 2 0/(2 cos 2 (j)—v 2 ) and the displaced phase satisfies tan 0(£) = — ^ tanh( 
where a = z>>/2 — v 2 and a = jr^a. The corresponding solution of the NLS equa- 
tion ([3]) with extremum at (£, r) = (0, 0) can then be written after some manipula- 
tions as: 



cosh(a^) — y 1 — 7ji> 2 cos(^r) 

This solution was found by Akhmediev et. al. (1987). They suggested a 
method of obtaining exact solutions of the NLS equation which is based on the 
substitution connecting the real and imaginary parts of the solution by a linear 
relationship with coefficients depending only on time. The method consists in 
constructing a certain system of ordinary differential equations, the solutions of 
which determine the solutions of the NLS equation. However, the explicit expression 
of the SFB already appeared in the previous papers (Akhmediev et. al. 1985; 
Akhmediev and Korneev, 1986). The SFB is also derived by Ablowitz and Herbst 
(1990) using Hirota's method (Hirota, 1976) and for v = 1 by Osborne et. al. 
(2000) using the inverse scattering technique. The SFB is the single homoclinic 
orbit to the cw as a fixed point of the NLS equation (Ablowitz and Herbst, 1990; Li 
and McLaughlin, 1994, 1997; Calini and Schober, 2002). The asymptotic behaviour 
for £ — > —00 of the SFB describes the linear instability of the cw solution for side- 
band perturbation (Akhmediev and Korneev, 1986). This modulational instability 
is known as the Benjamin- Feir instability in the context of water waves (Benjamin 
and Feir, 1967). A plot of the absolute value of this complex amplitude solution 
for 7*0 = 1 and v = 1/2 is given in Figure [TJa). 

2.2 Ma breather 

Secondly, we take £(r) = cosh(/ir), where \i = roy^/i. The differential 
equation for G becomes: 

d 2 T G = M 2 G - 3 M 2 ^G 2 - 2 M 2 ^-^G 3 . 

P P Z 

Comparing again with (TsJ) , the solution ^ is obtained with P{(j)) = — p 2 Q ((/))/ cos 0, 
Q 2 ((/)) = 2 cos 2 0/(2 cos 0+/i 2 ) and the displaced phase satisfies tan </>(£) = — tan(p£ 

where p = jl\/2 + p 2 and p = jr^p. The corresponding solution of the NLS (JiJ 
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(a) (b) (c) 



Fi gure 1: Plots of the absolute value of the SFB breather for v = 1/2 (left), the Ma breather 
for fi = 0.4713 (center) and the rational breather (right), for ro = 1 in all cases. For illustration 
purposes, the axis are scaled corresponding to /3 = 1 = 7. 



can then be written after some manipulations as: 

A(e,r) = A (0 ■ ( Z^i ^±M^ ML_ _ x \ . 

ycos(p^) — wlH- \ jl 2 cosh(/ir) J 

We call this solution Ma solution or Ma breather since it has been found by Ma 
(1979) using the inverse scattering technique of (Gardner et. al., 1967). The same 
technique is also used by Osborne et. al. (2000) to arrive at this solution for 
jl = y/2. A plot of the absolute value of the Ma breather for ro = 1 and jl = 0.4713 
is given in Figure [TJb). This choice of this value of jl is just for our convenience, 
which we will use again for the corresponding physical wave field. 

2.3 Rational breather 

Thirdly, we take ((r) = 1 — \v 2 t 2 . Similarly, the same solution can also 
be obtained by substituting £(r) = 1 + ^/i 2 r 2 . As a consequence, the differential 
equation for G now becomes: 

Comparing again with rt8l, we have P(0) = z> 2 / cos ^ QW = 1 + P 2 /( 4z>2 ) and 
the displaced phase satisfies tan (/>(£) = —2jrQ^. Substituting into the Ansatz Q, 
we obtain 

This rational solution is clearly different from the two solutions above, is derived 
by Peregrine (1983) as a limiting case from the Ma breather. He took the analytic 
expression and performed a double Taylor series expansion about the amplitude 
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peak (which occur at £ = = r) to arrive at his solution. Since this solution is 
confined in both space and time, some authors (Henderson et. al., 1999) also called 
it the 'isolated Ma soliton' and other authors (Akhmediev and Ankiewicz, 1997) 
called it the 'rational solution'. Because of its soliton-like feature in £, (Nakamura 
and Hirota, 1985) called it as an 'explode-decay solitary wave'. In this letter, we 
call it the 'rational breather'. A plot of the absolute value of the rational breather 
for ro = 1 is given in Figure [TJc). 



3. RELATION BETWEEN THE NLS BREATHER SOLUTIONS 

As shown by Dysthe and Trulsen (1999), there is a relationship among the 
breather solutions of the NLS equation. Different parameters are used for explicit 
expressions of the SFB and the Ma breather in Dysthe and Trulsen (1999) as well 
as in Grimshaw et. al. (2001). Introducing new parameters cp and $ G M by the 
following relations: v = V^sin^, a = sin2(/?, jl = V^sinh^ and p = sinh2^, we 
can write our expressions similar to the ones in these papers. 





v = m 


single soliton 






— > OO 




SFB 


Ma breather 





v -> 



Figure 2: The schematic diagram for the derivation of the Ma breather and the rational breather 
from the SFB. 



The SFB becomes the Ma breather if we substitute v = i/i and it becomes 
the rational breather if ^ — ^ 0. Similarly, the Ma breather becomes the rational 
breather for /i 0. Interestingly, the Ma breather becomes the one- soliton solution 
for ji — y oo (Akhmediev and Ankiewicz, 1997). Figure [2] explains the schematic 
diagram of these derivations. These relations can also be seen from the expression 
in the previous section. For the SFB, we take £(r) = cos(^r), which leads to the 
Ma breather if we substitute v = z/i, so that £(r) = cos(ifir) = cosh(/ir). Taking 
either v or \i — >> 0, we get ((r) = 1 — \ v 2 t 2 or £(r) = 1 + \h 2 t 2 , and both will lead 
to the rational breather. 

The relations have also consequence for the amplitude amplification factor 
(AAF), defined as the ratio between the maximum amplitude and the value of its 
background. The expressions for the breather solutions as given in this letter were 
shifted such that the maximum amplitude is at (£, r) = (0,0) and the value of 
the background is ro- For the SFB, the amplification is given by AAF s (z>) = 1 + 
V4-2z> 2 . For < v < the amplification is bounded and are 1 < AAF s (z/) < 3. 
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Figure 3: Plot of the amplitude amplification factor for the SFB and the Ma breather. The 
meeting point of the two curves is the AAF of the rational breather. 



The AAF for the Ma breather is given by AAF Ma (/i) = 1 + ^4 + 2jl 2 . Hence, for 
jl > 0, we have AAF Ma > 3. The AAF for the rational breather is exactly 3, which 
follows by letting the either v or fl go to zero 

AAF Ra = lim AAF S = 3 = lim AAF Ma . (9) 

In Onorato et. al. (2001), the AAF for the SFB is given as a function of the wave 
steepness and the number of waves under the modulation. The plot of the AAF 
for all the three breather solutions is given in Figure [3] For Jl — >• oo in the Ma 
breather, the single-soliton solution is obtained (Akhmediev and Ankiewicz, 1997). 

4. EVOLUTION IN THE ARGAND DIAGRAM 

In this section we present the evolution of the breather solutions in a complex 
plane. The real and the imaginary axes of this plane are denoted by the real 
and the imaginary part of the breather solutions after removing the plane-wave 
contribution, respectively. In this paper we refer this plane as the Argand diagram. 
Since the breather solutions depend on the space variable £ and the time variable 
r, the evolution curves are presented by parameterizing one variable and plotting 
for the other at different values. We will observe that the point (—1,0) plays a 
significant role with respect to all evolution curves, whether they are parameterized 
in space or in time. A study of this evolution curves for understanding experimental 
results on freak wave generation using the SFB solution has been presented recently 
(Karjanto, 2006). 

Figure ^ shows the evolution curves of the SFB for v = 1/ y/2. The parame- 
terized plot in position for different times is given by Figure |4^a). For this particular 
value of modulation frequency z>, all curves start from a point in the second quad- 
rant, which corresponds to £ = — oo, and end to a point in the third quadrant, 
which corresponds to £ = oo. During one modulation period < r < 27r/z>, we ob- 
serve that the curve reaches its maximum value for r = when the lm(F) vanishes 
and it reaches a minimum value for r = tt/u. The dotted circle indicates the plane 
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wave solution that has been removed. When we allow the modulation frequency 
v — >• 0, both the initial and the final points will coincide at (—1,0), as shown by the 
rational breather case in Figure^ a). The parameterized plot in time for different 
positions is given by Figure [4jb). All curves are straight lines and they are passed 
twice during one modulation period. For £ — >• ±oo, the lines shrink to a point in 
the third and the second quadrants, respectively. The longest line is reached at 
£ = when it lies on the real axis. Notice that ( — 1,0) becomes a center for all the 
lines. Moreover, for v — >> 0, the left-hand side of all the lines will meet at (—1,0), 
as shown by the rational breather case in Figure [6l 




i = 3tt 


/(2p) 


1 ^ 

5 = 




/(2/5) 





-2 2 4 

(b) 



Figure 5: The evolution of the Ma breather in the Argand diagram for £i = 2. 



Figure [5] shows the evolution curves of the Ma breather for ft = 2. The 
parameterized plot in position for different times is given by Figure |5|a). All the 
curves are elliptical, which indicate that they are ^-periodic. The largest curve 
occurs at r = 0, which has a circular form. As the time moves forward and 
backward, the ellipses are getting smaller until eventually shrink to a point at 
( — 1,0) for r = ±oo. By letting fl —> 0, the left-hand side of the ellipses at the real 
axis moves accordingly toward (—1,0), and eventually the pattern becomes the one 
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(a) (b) 
Fi glire 6 1 The evolution of the rational breather in the Argand diagram. 



as shown by the rational breather in Figure |6[a). The parameterized plot in time 
for different positions is given by Figure [5^b). Similar to the previous case, all the 
curves are straight lines and we notice clearly that they are centered at (—1,0). 
Different from the SFB solution case where the lines are at the right-hand side of 
( — 1,0), for the Ma breather case the point radiates lines to all different direction 
depending on the position. For one periodic position < £ < 27r/p, the longest 
line occurs at £ = and the shortest one occurs at £ = ir/p, being both of them 
lay on the real axis. When jl 0, the lines which are on the left-hand side of 
( — 1,0) will get shorter until eventually become the rational breather case as shown 
in Figure ^b). Figure [6] shows the evolution curves of the rational breather, which 
show the limiting cases for both v — ^ and jl —> 0. 



5. PHYSICAL WAVE FIELDS 

In this section, we show applications in surface water waves. Since the NLS 
equation is an envelope equation, the corresponding physical wave fields rj(x,t) are 
wave packets. Considering only first order contributions, the physical wave field is 
given by 

rj(x, t) = A(£ ry^ *-^ +cc (10) 

where c c denotes complex conjugate of the preceding term. The complex amplitude 
A is described in a moving frame of reference with £ = x and r = t — x/Vq. The 
wavenumber ko and the frequency ujq are related by the linear dispersion relation. 
For surface water waves, the normalized form of the dispersion relation is given by 
uj = Q(k) = kyj (tanh fc)/fc, the group velocity Vb = Q f (ko). 

Density plots for three breather solutions are given in Figure [7] For illustra- 
tions, we choose rn = 1, ko = 2ty (uo « y/2n) in all cases, modulation frequency 
v = 1/2 for the SFB and jl = 0.4713 for the Ma breather. For better representa- 
tion, the physical wave fields are shown in a moving frame of reference with suitable 
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chosen group velocity. The SFB has extreme values at x = £ = and is periodic 
in time. The Ma breather has extreme values at r = and is periodic in space. 
The rational breather is neither periodic in time nor in space, but isolated with 
its maximum at (x,t) = (0,0). Despite some differences, all the breather solu- 
tions show 'wavefront dislocation', when splitting or merging of waves occurs (Nye 
and Berry, 1974). A necessary condition for this phenomenon to occur is that the 
Chu-Mei quotient of the nonlinear dispersion relation at the vanishing amplitude is 
unbounded. When it occurs, the real- valued phase of the corresponding wave field 
is undefined at the vanishing amplitude, known as 'phase singularity'. More infor- 
mation on wavefront dislocation in surface water waves can be found in (Tanaka, 
1995; Karjanto and van Groesen, 2007). 




(a) (b) (c) 



Figure 7: Density plot of the physical wave field of the SFB for i> = 1/2 (left), the Ma breather 
for £i = 0.4713 (center) and the rational breather (right) corresponding to Figure [l] for ko = 2tt. 
The plots are shown in a moving frame of reference with suitable chosen velocity. 



6. CONCLUSIONS 

In this letter, we have shown that a description with displaced phase-amplitude 
variables, for which the phase is time independent, leads to the three breather so- 
lutions of the NLS equation: the Soliton on Finite Background, the Ma breather, 
and the rational breather. These are obtained as explicit solutions of the equa- 
tion for the displaced amplitude of the NLS equation. We described the relation 
between these three breather solutions and presented the amplitude amplification 
factor of each solution as function of the parameters. For the corresponding physi- 
cal wave fields, wavefront dislocation and phase singularity at vanishing amplitude 
are observed in all three cases. 
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